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Abstract: In new-physics processes that produce 6 or c jets, a measurement of the initial 
h or c-quark polarization could provide crucial information about the structure of the new 
physics. In the heavy-quark limit, the h and c-quark polarizations are preserved in the 
lightest baryons they hadronize into, and Ac, respectively. We revisit the prediction for 
the polarization retention after the hadronization process and extend it to the case of trans¬ 
verse polarization. We show how ATLAS and CMS can measure the 6-quark polarization 
using semileptonic A;, decays, and the c-quark polarization using A+ —)• pK~7r~^ decays. 
For calibrating both measurements we suggest to use tt samples in which these polariza¬ 
tions can be measured with precision of order 10% using 100 fb“^ of data in Run 2 of the 
LHC. Measurements of the transverse polarization in QCD events at ATLAS, CMS and 
LHCb are motivated as well. The proposed measurements give access to nonperturbative 
QCD parameters relevant to the dynamics of the hadronization process. 
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1 Introduction 

In order to fully explore the nature of new particles, both the sizes and the Lorentz struc¬ 
tures of their couplings will need to be measured. Probing the Lorentz structure is partic¬ 
ularly challenging as it often requires measuring the polarizations of final-state particles. 
Information about the polarization of colored decay products is typically washed away 
by hadronization. A well-known exception is the top quark [1, 2], which decays before 
it hadronizes. In this paper we show that, while challenging, the polarization of 5 and c 
quarks can also be measured at the LHC, despite hadronization. 

Knowing how to extract the 5-quark polarization could facilitate a variety of interesting 
measurements. For instance, in 5, —t- 55 decays one could examine whether the Higgs 
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coupling to b quarks has a CP-violating component, hlry^h, in analogy to the h —?• 
case [3]. Similarly, if a stop or a sbottom is discovered and its decay produces 6’s, one 
could determine whether it is the left-handed or the right-handed one, or, more generally, 
determine the left-right mixing angle. Also c quarks play an important role in a variety of 
new-physics scenarios, e.g. refs. [4-8]. 

As a proxy for the 6-quark polarization we are proposing to use the polarization. 
The Aft is a spin-1/2 baryon, which is produced in 6-quark hadronization both directly and 
from the decays of Eft and baryons, in comparable amounts. The main point is that, 
in contrast to the B mesons, the Aft is expected to retain the polarization of the 6 quark 
to a high degree, at least in the heavy-quark limit [9-11]. About one out of ten 6 quarks 
produces a Aft, and these events can be used for extracting the 6-quark polarization. 

We define the fraction of polarization retained in hadronization to a Aft as 


_ ^(Ab) 

'P~ V{b) ’ 


( 1 . 1 ) 


where V{b) is the polarization of the 6 quark as it exits the hard process and V{Ab) is the 
Aft polarization when it decays. In general, rp depends on the initial polarization direction, 
V{b). If the 6 is either longitudinally or transversely polarized, then rp is a number, or 
ry, respectively, while it is a tensor in general. In the heavy-quark limit, rrn, ^ Aqcdj one 
has rp = 1. For the physical 6 mass we thus expect to be 0(1) [9-11], where the precise 
number depends on relatively uncertain hadronization parameters. We suggest to measure 
rp at the LHC in Standard Model (SM) processes with polarized 6 quarks. The results 
will allow interpreting similar future measurements of 6-quark polarization in new-physics 
processes. 

As long as the hard scale, Q, at which the 6 quarks are produced is much larger 
than the QCD scale, Q ^ Aqcd, tbe 6-quark hadronization and the subsequent evolution 
factorize from the short-distance production process. Therefore rp is a universal quantity, 
independent of the exact mechanism that produces the initial 6 quark. In general, rp 
depends on the scale, Q, and the fraction of the 6-quark momentum carried by the Aft, z. 
The important point is that once we know rp{z) at a given scale, we can calculate it at 
a different scale using the known renormalization group (RG) evolution of fragmentation 
functions. A measurement of rp{z) using a SM process at some scale Q will then enable us 
to know rp{z) at any scale and use it in new-physics measurements. Moreover, the effects 
of scale dependence are small as long as the characteristic scales of the measurements are 
similar. Thus, at the first stage, measurements inclusive in z are sufficient. Only once we 
enter a precision era will one need to take into account the effects of running. 

Depolarization occurs both during and after hadronization. During hadronization the 
flip of the 6-quark spin occurs via QCD-scale processes. It is Aqcd/ mb suppressed because 
the 6-quark chromomagnetic moment is ^b oc l/mft and is, as such, small. After hadroniza- 
tion, depolarization occurs mainly because Aft’s are also produced from ^ decays whose 
lifetimes are longer than the timescale for hadronization into distinct mass eigenstates Eft 
and E^, i.e., T (*) < ms* — w-S;,- Even though this effect vanishes in the formal nib —t oo 

2 j , b 
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limit, it is 0(1) for the physical 6-quark mass [11]. The dominant depolarization effect is 
therefore due to the ^ decays. 

Evidence for longitudinal polarization in Z —)• 66 decays has already been seen at 
LEP [12-14], but precise measurements of were impossible. At the LHC, the Z —)• 
bb sample suffers from a large QCD background, pp —)• 66 -|- X [15], which makes the 
measurement difficult despite the fact that the background 6’s are only slightly (and just 
transversely) polarized. In contrast, as we demonstrate in this paper, the 6’s from top- 
quark decays at the LHC allow for a clean measurement of at ATLAS and CMS with 
the upcoming Run 2 datasets. 

It would also be useful to measure using the transverse polarization of 6’s produced 
in QCD events. The polarization in QCD events arises at NLO and for large momenta 
behaves like V{b) ~ agnib/pb, where pb is the 6-quark momentum [16]. Since it is larger 
for softer 6 quarks, the corresponding A;, decays are probably easiest to reconstruct at 
LHCb (although to use A;, as a 6-quark proxy, the 6 quarks still need to be hard enough for 
factorization to apply). However, the polarization varies significantly as a function of the 
parton-level kinematics of the event, and even changes its sign for some of the contributing 
processes [16]. The limited angular coverage of LHCb may hinder using this kinematic 
dependence, which is ignored in the existing LHCb measurement [17]. Therefore, low- 
Pt measurements by ATLAS and CMS, e.g. along the lines of refs. [18, 19], seem to be 
motivated as well. 

An additional motivation for measuring the A;, polarization (and a few related quanti¬ 
ties, as we will discuss) in SM processes is that it can teach us a lot about the hadronization 
process and provide access to several nonperturbative QCD parameters. As we will review, 
the present knowledge of the relevant physics is incomplete. The results of the measure¬ 
ments can also be useful in tuning Monte Carlo generators. 

In the case of c quarks, the physics of the relevant baryons (Ac, Sc*^) is qualitatively 
similar to the 6-quark case. It is likely that an 0(1) fraction of the polarization is preserved 
despite the fact that rric 3> Aqcd is not a very good assumption. The transverse polariza¬ 
tion of Ac’s from QCD production has already been seen in the fixed-target experiments 
NA32 [20] and E791 [21], but theoretical interpretation is difficult because soft QCD effects 
may play a major role for the relatively low Ac momenta probed in these experiments. We 
will discuss how can be measured at ATLAS and CMS using a ti sample, in which 
polarized c quarks are available from W~^ —)• cs decays. 

The rest of the paper is organized as follows. In section 2 we describe the basic 
properties of the baryons of interest, while in section 3 (and appendices A and B) we 
study the polarization transfer from the heavy quark to the baryon. In sections 4 and 5 
we analyze how to measure the polarization of the relevant 6 or c baryons at the LHC 
and propose specific analyses for such measurements in pp —)• ti. In section 6 we discuss 
how to obtain additional information by studying , Sc contributions in isolation. We 
summarize in section 7. Appendix C describes the relation between rp and fragmentation 
functions. 
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2 Bottom and charmed baryons 

A b quark can combine with a light diquark^ to produce a baryon. Most commonly, the 
diquark is made out of u and/or d quarks, producing either the isosinglet spin-1/2 baryon, 
Afc, or one of the isotriplet spin-1/2 or spin-3/2 baryons, S;, and respectively. The latter 
decay primarily through QCD as A^ tt, while —)• Tij, decays can be neglected. 

The Ajj decays via weak interactions and can be treated as an asymptotic state in our 
discussion. 

The probability for a b quark to fragment into any baryon is 

/baryon = (8.0 ± 1.0) % , (2.1) 

based on LEP measurements of Z ^ bb decays as summarized in table 5 of ref. [23]. This 
number includes /baryon = /aj, + fst + where and Qb are baryons that contain one 
and two strange quarks, respectively. Baryons that decay to A;, before the b itself decays, 
such as are included in (and similarly for f^b and /o;,)- We estimate the Ab 

contribution to /baryon to be about 85% [24, 25], while the rest is primarily Eb, which is 
studied in appendix B. This estimate is obtained from the relative rates of the 6 —)• A;, —)• 
J/ipA, b —)• E^ —>■ J/'tJjE~, b —>■ processes measured in QCD events at 

the Tevatron [24], using theoretical predictions for the branching ratios to J/ijj [25], and 
assuming /■=o = /—. For numerical estimates in the rest of the paper we will therefore use 

-f, 

/a. = 7% . (2.2) 

In the near future, the LHC experiments will likely shed more light on the baryon frag¬ 
mentation fractions. 

The c quark has a similar spectrum of baryon states. The fragmentation fraction of a 
c quark into a Ac, based on LEP measurements [26], is 

/a, = (5.7 ±0.7)%. (2.3) 

Several experiments reported that in QCD events /a^ varies significantly as a function 
of the 6-hadron px, even for px Aqcd, where factorization is expected to hold [23, 27- 
30]. This should not be interpreted as a variation of /a^ from eq. (2.2) with the energy 
scale of the process. Events with the same px of Ab can come from b jets with very different 
px of the original b quark, by which we mean the total px of the b jet, after adding the 
reconstructed neutrino px if relevant. One gets contributions from b jets where the A^ 
carries most of the 6-quark momentum as well as from much harder 6 jets where the Ab 
carries only part of the momentum. Because the QCD production cross section changes 
rapidly with the 6-quark px, a small difference in the shapes of the fragmentation functions 
of different 6 hadrons can translate into a large difference in their contributions to fixed 
hadron-p-r bins (see also ref. [31]). This can lead to an apparent px dependence of the 
fragmentation fractions even if it is absent at the fundamental level. As discussed in more 

^The concept of a diquark, as the state of the light degrees of freedom within a heavy baryon, has precise 
meaning in the framework of the Heavy Quark Effective Theory (HQET). For a review, see, e.g., ref. [22]. 
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detail in appendix C, a much clearer interpretation would be obtained if the measurements 
were performed in terms of fixed reconstructed 6-quark px rather than 6-hadron pT- In that 
case one expects to see only a slow (logarithmic) dependence on the hard scale due to RG 
evolution. That is, it would be desirable for the LHC experiments to perform measurements 
of the differential cross sections da/dpT in terms of the reconstructed 6-quark pT- When 
enough data are available one should also perform measurements of d?a/dpxdz, where z is 
the Af, momentum fraction relative to the initial pT of the 6 quark. 

As mentioned above, the A;, polarization carries information about the initial 6-quark 
polarization and the leading depolarization effects are due to S;, and decays. To describe 
the relative production probabilities of A^, S;,, and we write their wave functions in 
terms of diquark and 6-quark states. The diquark can be a spin singlet, S', or a spin triplet, 
T. This allows for four possible spin conhgurations. So ,T+i,ro,r_i, where the subscripts 
denote the spin projection along the spin-quantization axis. Using the same quantization 
axis for the spin of the 6 quark, the baryon mass eigenstates are 

\\±l) = \b^)\So), (2.4) 

1^6,±l) = ± yi \b^l)\T±i ), (2.5) 

= \ll Kl)\To) + |6^i)|T±i), = |6±i)|r±i) . (2.6) 

The relative probabilities to produce So and To,±i around the 6 quark control the 
relative size of direct A;, production and its production from decays of various states. 
These probabilities can be parameterized in terms of two nonperturbative parameters, 
0 < A < oo and 0 < rci < 1 [11], 

P[T_i] and P[r+i] are equal because QCD is parity invariant. The parameters A and wi 
are inclusive over the momentum fraction z of the A;, inside the 6 jet. They do, however, 
have a weak dependence on the hard scale, Q, as discussed in appendix C. In the remainder 
of this section we discuss what is known about the values of A and wi. 

The parameter A is the ratio of the production rate and the direct A;, production 
rate. While the CDF collaboration has measured the masses and widths of the [32, 33], 
it has not determined their production rates. We therefore estimate A using the statistical 
hadronization model (for a brief overview, see ref. [34] ), according to which the production 
rate per degree of freedom is proportional to 

, ( 2 . 8 ) 

where m is the mass of the hadron and T ~ 165 MeV [34]. This gives 

A~2.6, (2.9) 

for both the bottom and the charm systems. 
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The value in eq. (2.9) is significantly larger than the estimate in ref. [11], which set 
A = 9 PARJ(4), where PARJ(4) is the Pythia6 parameter in the Lund fragmentation 
model describing the probability for forming a spin-1 vs. spin-0 diquark [35, 36], and the 
factor of 9 is the multiplicity ratio of isotriplet spin-1 and isosinglet spin-0 diquark states. 
The equivalent PythiaS parameter is StringFlav:probQQ1toQQ0 [37]. Depending on 
the choice of the Pythia tune [36, 37] this gives values of A between 0.24 and 0.45. The 
discrepancy with the estimate in the statistical model is likely due to the fact that the 
Pythia tunes are based on light hadrons. There is no reason to expect this phenomeno¬ 
logical parameter to have the same value for heavy-quark hadrons. On the other hand, 
the quark-diquark model of heavy-baryon production in ref. [38] predicts A ~ 6 for both 
the bottom and the charm systems. Though somewhat larger, this is of the same order 
of magnitude as our generic estimate in eq. (2.9). The measurement of the relative Sc/Ac 
yield by E791 [39] gives a somewhat smaller value than eq. (2.9), A ~ 1.1 (in extrapolating 
to S* we included the factor R from eq. (3.31), discussed below). The measurement of Si, 
and production by DELPHI [40], in combination with eq. (2.2), gives 1 < A < 10, again 
favoring eq. (2.9) over the Pythia parameter. 

The parameter wi accounts for the possibility that the fragmentation axis breaks the 
rotational symmetry in the spin-1 diquark production. The isotropic case is when wi = 2l‘i. 
DELPHI studied the angular distribution of ^ Afovr decays at LEP [40-42] finding 

=-0.36 ±0.30 ±0.30. (2.10) 

Since negative values of wi are not physically meaningful this suggests that rci ~ 0. In 
contrast, an analogous measurement in the charm system by CLEO at CESR gave [43] 

rci = 0.71 ±0.13, (2.11) 

consistent with the isotropic case. A theoretical calculation [38] based on a quark-diquark 
model gives wi ~ 0.41 and wi ~ 0.39 for the bottom and charm system, respectively. 
The uncertainties on these estimates due to assumptions made in ref. [38] may be large. 
Eor instance, finite-width effects, describing the interference between Si, and S^, are quite 
important (cf. section 3.2.2), but were neglected in ref. [38]. Eor these reasons, we shall 
treat wi as a yet-unknown parameter. Eor other discussions of rci, and an analogous 
parameter W 3/2 relevant to excited mesons, see refs. [11, 44-49]. 

3 Af, polarization and ^ decays 

When the b quark emerges from the hard process, it loses only about 2asl‘i'K ~ 3% of 
its polarization to gluon radiation [50]. During the fragmentation process, in the exact 
heavy-quark limit, ruft/AgcD cOj the QCD interactions cannot change the spin of the b 
quark because its chromomagnetic moment is proportional to Ijmi,. This is the case for 
all b hadrons. The additional special property of the A^ is that in the heavy-quark limit 
its light degrees of freedom form a spin-0 state, and thus do not affect the spin of the b 
throughout the A^ lifetime. 
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Parameter 

(MeV) 

Parameter 

(MeV) 

m-Si, - ruAi, 

194 ±2 

ms, - mA, 

167.4 ±0.1 

ms* - ruAi, 

214 ±2 

ms* - mA, 

231.9 ±0.4 

A = ms* - mSi, 

21 ± 2 

A = ms* - ms. 

64.5 ±0.5 

Ts, 

7±3 

Ts, 

2.2 ±0.2 

Ts* 

9±2 


15 ± 1 


Table 1. Measured charge-averaged masses and widths of (left) and (right) [51]. is 
also known as Sc(2455), and S* as Sc(2520). 


As pointed out in ref. [11], the dominant depolarization effect is that the final A;, 
sample contains contributions from 6’s hadronizing into ' that subsequently decay to 
Aft. In the , depolarizing chromomagnetic interaction between the spins of the b quark 
and the diquark acts over relatively long timescales given by the ' lifetimes. We have 
r (*) < A < Aqcd, where 

A = ms* - msi, (3.1) 

is the hyperfine splitting, see table 1 (left). Therefore, hadronization to distinct mass 
eigenstates occurs before they decay. Since some of the states are not eigenstates 
of the 6-quark spin, see eqs. (2.5)-(2.6), the depolarization effect can be of 0(1). 

This effect vanishes in the mb —?• oo limit. In this limit, the decay widths T (*), given by 
the HQET expression in eq. (3.19) below, remain largely unchanged because m (*) —m^^ is 
approximately independent of mb- The hyperfine mass splitting, on the other hand, scales 
as A oc Ijmb, so that for large enough mb one has T (*) 3> A and no depolarization occurs. 
However, this is not the situation realized in nature. 

In the rest of this section we describe the ' production and decays and how these 
influence the Aft polarization. We show that the polarization of Aft’s from decays 
depends on both the magnitude and the direction of the original 6-quark polarization. The 
results will be expressed in terms of the angle 6 p, defined in the ' rest frame as the 
angle between the initial 6-quark polarization and the fragmentation axis, which lies along 
the direction of motion of the 6 quark (see figure 1). For 6 quarks from top or Z decays, 
the electroweak interaction produces longitudinal polarization, namely 9p = 0. This was 
the only case analyzed in ref. [11]. For 6 quarks from QCD production, where a small 
polarization arises at NLO [16], Op = 7r/2. In new-physics models. Op can in principle 
have any value. For instance, 6 quarks produced in decays of a right-handed sbottom to a 
bino will have a longitudinal polarization of -|-1. Transversely polarized 6’s can arise, for 
example, due to a broad resonance interfering with QCD processes, similar to what has 
been discussed in the context of the transverse polarization of top quarks in ref. [52]. 

The Ac-Sc system is described by qualitatively the same physics as the Aft-S^ ^ 
system. The parameters of the relevant baryons are shown in table 1 (right) and the 
corresponding results for the polarization will be presented in section 3.4. 
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fragmentation axis 
(direction of motion) 


1 r{b) 



Figure 1. The angle Op and the polarization retention factors and tt- 


(*) 

3.1 Production of SI ^ and their decays 


We orient our coordinate system such that the b polarization axis in the S^ ' rest frame 
is pointing along the z axis. The parameterization of production probabilities in eq. (2.7) 
applies to the spin states of the spin-1 diquark along the fragmentation axis, These 

are expressed in terms of the states along the b polarization axis, \Tm), as 


i^m'(^p)) = '^Rm'm{Gp) \Tm) , 


(3.2) 


where 




/ 2 ^P ■ 2 ^P \ 

2 V2 2 

sin 6 p sin 9p 

cos Op — - 


sm 


%/2 

2 


V2 


Up 2 '^p 

COS 


(3.3) 


2 V2 2 / 

for m,m' = —1,0,-|-1. Combining a 6-quark state with spin along the z axis with a 


diquark spin state we obtain 

\b^i)\r^,{0p)) = Y,Rm'.^{.ep) 


M 


M 


(3.4) 


(*) 

where M is the spin component of the ' along the 2 ; axis 
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(*1 

In the heavy 5-quark limit, the decays —)• A^tt proceed via the decay of the internal 

spin-1 diquark, Tq^±i, to the spin-0 diquark, Sq, leaving the b quark and its spin unaffected. 
Since the initial diquark has spin 1, while the final diquark and the pion are spinless, the 
orbital angular momentum state of the decay products must be £ = 1. Therefore, a ^ 
spin state described by J, M decays to a state of the form 

|^(J,M))oc /dcos^dc^ l,M-s| J,M)yi^-®(0,0)|0,(/.)|s). (3.5) 

S 

Here, 0, (j) describe the direction of motion of the pion in the ^ frame, s is the spin 
along the 2 : axis, and Y^{6^ (j)) are the spherical harmonics. 

3.2 Effect of ^ decays on Af, polarization 
3.2.1 Aft polarization in the limit of narrow 

For simplicity, we first assume that the ^ widths, r^(*), can be neglected relative to 
the mass splitting A = m^* — rriY,^. In this case and Sft decay incoherently since the 
different pion energies in their final states prevent interference. Taking into account the 
amplitudes for producing the various spin states based on eq. (3.4) and the decay 
amplitudes from eq. (3.5), an initial state \h_^i)\T!^i{6p)) produces the state 

I'k) OC / d cos 6 d(j) ^ 5) + 5! Ij ''TT-) X 

m M 

X s; 1, M - s I J, M) cj,) \e, ct>) \s) , (3.6) 

S 

with J = 5 and | for the Eft’s and E^’s, respectively. We shall assume that the pion degrees 
of freedom \0,4>) will not be used in the measurement due to experimental difficulties 
discussed in section 6. By tracing over the pion degrees of freedom we readily obtain the 
density matrix of the Aft spin 

oc Tre,^ |T) (T| , (3.7) 

where T = Eft or E^. 

The total density matrix, combining both Eft and E^ decays, is given by 

p oc E pq,p^ , (3.8) 

where is the probability to produce a particle of type 'k. From eq. (3.4), ps*/pSj, = 2. 
This factor receives a small correction from the fact that a heavier state is less likely to be 
produced in fragmentation — the Boltzmann factor in eq. (2.8) suppresses E^ production 
relative to Eft production by a factor of 

i? = ~ 0.88 . (3.9) 

The deviation of R from unity is an 0{AqcD/'mb) effect and we have been neglecting 
other effects that are formally of the same order. However, keeping R ^ 1 will facilitate 
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comparison with the results of the next section, where we go beyond the narrow-width 
approximation. Furthermore, measurements in the D-D* system, which is analogous to 
the system [ 11 ], point to the phenomenological relevance of i? 7 ^ 1. As discussed 

in ref. [53] and references therein, the well-measured deviation of the D/D* multiplicities 
ratio from the naive prediction is in agreement with the expectation from the statistical 
hadronization model. At the same time, the spin alignment of D* mesons is in agreement 
with expectations from Clebsch-Gordan coefficients without requiring l/ruc corrections [54]. 
Combining the production probabilities from eq. (3.4) with this additional correction factor, 
we rewrite the total density matrix as 


/O oc pEt + PS* ■ (3.10) 

As a last step we average the contributions to p from all diquark spin components m' 
with relative probabilities determined by the parameter wi from eq. (2.7). 

Finally, we normalize the density matrix to Trp = 1 and use the relation 

(3-11) 


to determine the polarization V. By symmetry, the polarization in our case can only lie 
in the xz plane, the plane formed by the initial h polarization and the fragmentation axis. 
Eq. (3.11) is thus explicitly 

P = 5 [(1 + P.) It) (tl + (1 - P.) It) (tl + P. (It) (tl + It) (tl)l . (3.12) 

The two components of the polarization vector are 

_2R- I + 2{l + R)wi + {I + R) {2 - 3wi) sin^ Op 

" ~ 3(l + 2i?) 

1 + i? / 2 \ 

^ 1 + 2R ( ~ 3 ) ■ 


(3.13) 

(3.14) 


Above we included only A^’s produced from T,/ ^ decays, while directly produced A^’s will 
be added below. 

For generic 9p the polarization vector changes direction relative to the polarization of 
the original b. This means that rp in eq. (1.1) is a tensor in general. However, if the 
initial 6 -quark polarization axis and the fragmentation axis are collinear, = 0 or vr, or 
are orthogonal to each other. Op = 7r/2, the polarization direction remains unchanged, as 
expected by symmetry. A longitudinally polarized 6 quark therefore results in a longitu¬ 
dinally polarized A;, and a transversely polarized 6 quark in a transversely polarized A;,. 
For isotropic diquark production, wi = 2/3, the magnitude of the final polarization is 
independent of Op and its direction is unchanged, as expected. 

For a longitudinally polarized 6 quark. Op = 0, the general result in eq. (3.13) reduces 


to 


ryL 

' Z 


2R-\^2{\^R)wx 
3(1 -h2i?) 


~ 0.09 -|- 0.45t(;i, 


(3.15) 
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and for a transversely polarized b quark, = 7r/2, to 


' Z 


4i? + 1 — (1 + i?) Wi 
3(1 + 2R) 


~ 0.55 - 0.23u)i. 


(3.16) 


Including the direct A;, production from fragmentation, the corresponding polarization 
retention factors from eq. ( 1 . 1 ) are 


rL,T 


l + AVz’'^ 
1 + A 


(3.17) 


3.2.2 Aft polarization for finite ^ widths 

The widths are only two to three times smaller than their mass splitting, cf. table 1. 
Sizeable interference effects may thus be present, so we extend our calculation to the case 
of finite widths. After the production of a Sft~S^ superposition state with energy E, and 
its decay to Aftvr, the state vector is 


1^;) oc f dcos9d(f) Rm'mjdp) M | 5 ,+ 5 ; l,m) 

m J^M 


Pn{E) 


E — mj + iT{E)/2 


X s; 1, M - s I J, M) cP) \9, cP) \s ). 


(3.18) 


Here mj is the mass of Eft or for J = ^, |, respectively. The pion-momentum factor 
p^{E) ~ ~ derives from the pion coupling in heavy-baryon chiral per¬ 

turbation theory [22, 55]. Correspondingly, for the width function r(ili) in the propagator 
we use 

= ( 3 - 19 ) 

where /jr — 93 MeV. This should satisfy r(m {*)) ~ T (*). We take the axial-vector current 
coupling gA to be 0.63 instead of 0.75 measured in neutron decay, to better reproduce the 
measured ' and Ec ' decay widths (see table 1). This choice corresponds to Tsj, ^ 
6.1 MeV, Ts* — 10.7 MeV, Ts^ — 2.1 MeV, and T^j — 15.7 MeV. 

We then proceed as in the previous section. From the density matrix 


p{E)^TTe,<p\E) {E\ 


(3.20) 


we find the polarization of Aft’s produced from E^ '’s with energy E to be 
V^{E) = 1 - 2(2 - wi)fiE ), rj iE) = l-{2 + wi)f{E ), 

in the longitudinal and transverse case, respectively, where 


fiE) = 


4 (ms* -msj^ 


12 [2{E - msJ2 + {E- ms*)^] + 9r2(^) ' 

The resulting behavior is shown in figure 2 for the case of longitudinal polarization 


(3.21) 


(3.22) 
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f*) f*") 

Figure 2. Polarization of Af,’s produced from ^ decays as a function of the ’ energy E. The 
polarization (red curves) is shown for the longitudinal case with wi = 0 (solid), 2/3 (dashed) and 
1 (dotted). The Ef, and production peaks are shown in black (arbitrary y scale). Vertical lines 
show ±r„(.) /2 ranges around the nominal masses, and horizontal lines indicate the values of the 
polarization in the narrow-width limit. 


Since we assume that the pion is not identified we average the polarization over E. 
The corresponding density matrix is 

poo 

poc dEP t^{E) p{E), (3.23) 

where the p-k{E) factor accounts for phase-space integration; this is in addition to two such 
factors already present in p{E) via eq. (3.18). The Boltzmann factor with T ~ 165 MeV 
is the equivalent of eq. (2.8). 

Numerically, we find the polarizations of A^’s from E/ decays in the longitudinal and 
transverse scenarios to be 

T’f'~ 0.23-b 0.38r(;i, -pj ~ 0.62 - 0.19u)i. (3.24) 

These results should be compared with eqs. (3.15) and (3.16) that were derived in the 
narrow-width approximation. We see that finite-width effects are non-negligible. 

The overall polarization retention factors from eq. (1.1), as computed from eq. (3.24) 
using eqs. (3.17) and (2.9), are 

ri ~ 0.45, 0.63, 0.72, rr ^ 0.72 , 0.63 , 0.58 , (3.25) 

for wi = 0, 2/3, 1, respectively. If we allow the parameter A to differ from our estimate 
in eq. (2.9) by up to a factor of two, the ranges of possible values of and vt extend to 
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(3.26) 


0.36 < < 0.78 and 0.52 < rr < 0.78. The polarization for arbitrary Op is given by 

Pz = vl cos^ Op + rT sin^ Op, Px = {tl — ^t) sin Op cos Op . 

In appendix A, we derive approximate analytic expressions that describe the results 
we obtained here. We also present an alternative picture of the physics, in which the 
depolarization happens due to oscillations between 6-spin eigenstates, analogous to K^~K 
oscillations. 


3.3 Results from LEP 


The Afe polarization has been measured, although with a large uncertainty, in Z decays at 
LEP, using the semileptonic decays of the A;,. The polarization of 6’s produced in Z decays 
is expected to be longitudinal and given by 

‘2VjyCll) 


-0.94, 


(3.27) 



where u;, = — 1 | sin^ 0 ^ and aj, = — 1 are factors in the vector and axial-vector couplings 

of the Z to 6 quarks. QCD corrections reduce this value by about 3% [50] . 

ALEPH and DELPHI used the variable {Ei)/{Ei,) proposed in ref. [56] (for a review 
of earlier literature on the subject, see ref. [57]), obtaining 


P(Ab) = -0.23tolt (stat.) (syst.) (ALEPH [12]), (3.28) 

P{Ab) = -0.49+|5Jo (stat.) ±0.17 (syst.) (DELPHI [14]), (3.29) 

while OPAL used a fit to the E^/Ey distribution, obtaining 

p(Afe) = -0.56 (stat.) ± 0.09 (syst.) (OPAL [13]). (3.30) 


Even though the precise value of the polarization retention factor cannot be determined 
from these results due to the large uncertainties, they do suggest that some polarization 
loss is present (i.e., r/, = 1 is excluded), but still is 0(1). Both facts are in agreement 
with expectations, see eq. (3.25). Large values of wi seem to be disfavored, especially by 
the ALEPH result. 


3.4 The charm case 

The ideas of this section can also be applied to c quarks. Similarly to eq. (3.9) we have 

Ra, ^ 0.68 , (3.31) 

with which we find the Ac polarizations from Sc decays in the longitudinal and transverse 
scenarios to be 

(pf)^^ ~ 0.07±0.46r(;i, ~ 0.54 - 0.23r(;i, (3.32) 

to be compared with eq. (3.24) for the A;,. The total polarization retention factors for 
wi = 0, 2/3, 1 are 

(ri)^^ ~ 0.33, 0.55, 0.66, (rr)^^ ^ 0.66 , 0.55 , 0.50, (3.33) 
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respectively. If we allow the parameter A to differ from our estimate in eq. (2.9) by 
up to a factor of 2, the ranges of possible values extend to 0.22 < ^ 0.73 and 

0.42 < (rr)^^ < 0.74. 

An important caveat is that 0 (Aqcd/?^c) corrections are likely to be larger than the 
0 (Aqcd/?^ 6) corrections that we have been neglecting in the b system. In particular, 
it may no longer be a good approximation to neglect the polarization loss in the initial 
stage of the fragmentation occurring at the QCD timescale. Nevertheless, even with these 
effects, the polarization retention factors are likely to remain 0(1). This is supported by 
the observation that even the A’s produced in Z —)■ jj decays at LEP retain an 0(1) 
fraction of the strange-quark polarization [58-60]. It should be noted that much of the 
polarization reduction in the case of A’s at LEP is not due to polarization loss during 
the s-quark hadronization, but because of an 0(1) contamination from unpolarized A’s 
produced from s quarks appearing in the fragmentation process [60]. Such contaminations 
are expected to be smaller in the Ac case. A large transverse Ac polarization was measured 
in QCD processes in the fixed-target experiments NA32 [20] and E791 [21], but theoretical 
interpretation of these results is difficult (see also ref. [61]) because the typical pr’s of the 
Ac’s (~ 1.5 GeV) were not much larger than the QCD scale. 

4 6-quark polarization measurement via semileptonic A;, decays 

Here and in the next section we outline several possible strategies for A;, and Ac polar¬ 
ization measurements in ATLAS and CMS. The ultimate goal is to study b- and c-quark 
polarizations in new-physics processes. As a SM calibration we propose the ti sample. 
The top decay acts as a “standard candle”, fixing the polarization retention factor of 
b quarks (from primary top decay) and of c quarks (from W decay). In both cases, the 
polarization of the initial quark is to a good approximation completely left-handed, i.e., 
V(b) ~ —1, V(c) ~ —1 in our convention. 

4.1 Properties of the decay 

To measure the Af, polarization one can use its inclusive semileptonic decay 

Ah^X.ru, (4.1) 

proceeding via the partonic b —?■ cW~* —>■ c£~v transition. Here, Xc is an inclusive final 
state with nonzero charm quantum number. The branching ratio is H(Af, —)• Xc(-~i>) ~ 10% 
for each lepton flavor [51]. The kinematic distributions of the charged lepton and neutrino 
in eq. (4.1) have been obtained using operator product expansion and HQET, and are 
under good theoretical control [62]. They are 

=\(1 +otiV (Ab) cos 6 i), i = ioicn, (4.2) 

1 Ab d cos bi 2 

where 9i (9y) is the angle in the A^ rest frame between the lepton (neutrino) momentum 
and the A;, polarization. The distribution is uniform in the azimuthal angle (pi ((pi,). At 
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leading order in and ag, the coefficients ai^i, multiplying the A;, polarization, 

sometimes called the spin-analyzing powers or the decay asymmetry parameters, are 

— i + 4xc + 12 x^ — + 12 x^ log Xc + Sx^ log Xc 

Oil = —^^—q-1-9- 

1 — 8xc + SXg — x^ — 12 Xc log Xc 
«!/ = !, (4.4) 

where Xc = There are no corrections to eqs. (4.3) and (4.4) at ©(Agco/^^fe), while 

©(Aqqj-j/tti^) corrections [62] are negligible for our purposes. Higher-order corrections in 
as are also small; they increase ai by ~ 5% and decrease a^, by ~ 1% [63, 64]. 

For longitudinally polarized b quarks, the angles 6 i and 9^ should be measured with 
respect to the A^ flight direction. This is the case for b quarks from Z and top decays 
and in many new-physics models. For b quarks from QCD production, the polarization is 
perpendicular to the plane formed by the b quark and colliding partons [16]. 

Afc’s of opposite polarization give the same distributions as eq. (4.2). This means that 
in Z or tt events, for example, the decay products are distributed in the same way relative 
to the 6 -jet axis regardless of whether the jet originates from an initial 6 or 6 quark. 

We note that the neutrino is more sensitive to the A^ polarization than the charged 
lepton, see eqs. (4.3), (4.4). The polarization measurement requires knowing the A;, rest 
frame and thus it is necessary to reconstruct the neutrino momentum regardless of whether 
it is used as a spin analyzer. Another benefit of using the neutrino is that inclusively ckj/ 
is very close to maximal. Therefore, it must remain close to 1 also if we restrict the 
analysis to a not-too-small subset of the semileptonic decays. This is advantageous since 
different semileptonic decay modes or kinematic regions may have different efficiencies, 
either due to experimental limitations or due to cuts applied for background reduction. 
An important intrinsic background arises from semileptonic H-meson decays. Even though 
these decays are isotropic, their presence in the sample dilutes the observables sensitive to 
the Af, polarization. 


~-0.26, (4.3) 


4.2 Strategy for Aft-polarization measurement 


We suggest to measure the forward-backward asymmetry of the neutrino, Afb, in the Aft 
rest frame along the expected direction of polarization. 


Afb 


N+ - N_ 
N+ + N_ ' 


(4.5) 


Here, A^+ and N- are the numbers of events with cosOy > 0 and cosO^, < 0, respec¬ 
tively. The neutrinos in signal events are distributed according to eq. (4.2) while in the 
semileptonic decays of B mesons they are distributed isotropically. As long as the H-meson 
decays are reconstructed correctly they simply dilute the asymmetry. Afb then measures 
the polarization as 

P(A.) = ^, 

where / is the signal event fraction. The statistical uncertainty on Afb is 



AAfb 




(4.7) 
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Figure 3. An example tt event that can be used for measuring the polarization of b quarks produced 
in top decays. 


where N = + N- is the total number of events. 

In the rest of this subsection we propose how to tag b jets with Af, —)■ Xcii', recon¬ 
struct the Xc, and reconstruct the neutrino. Here, we keep the discussion general, but in 
subsection 4.3 we will analyze, as an explicit example, the polarization measurement of b 
quarks produced in top decays, illustrated in figure 3. We will estimate the sensitivity of 
the proposed strategy using efficiencies of similar procedures available in the experimental 
literature. 

We focus on measurements of 'P(Ab) inclusive over the A;, momentum fraction, This 
is sufficient as an initial calibration and is the only type of measurement that needs to be 
performed on new-physics samples. The next experimental step would be measuring V{Ab) 
in SM calibration samples in bins of z, which will provide inputs to the RG running of the 
polarization retention factors as explained in the introduction and appendix C. 

4.2.1 “Soft muon” b tagging 

Most LHC analyses apply standard 6-tagging algorithms based on the lifetimes of the 6- 
ffavored hadrons and/or the 6-quark mass. A better choice for our purposes is a “soft muon” 
6 -tagging algorithm. It demands a muon among the jet constituents, where the muon need 
not be isolated (unlike the hard lepton from the t —)• IT —)• f chain in section 4.3 below). 
The muon impact parameter and its transverse momentum with respect to the jet axis, 

give additional discrimination from non-6-flavored jets. 

The reason why “soft lepton” 6 tagging is not a popular choice in high-p-p analyses 
is that its efficiency is limited by the sum of the b ^ i and b ^ c ^ i branching ratios 
(roughly 10% each [51], times two lepton flavors). For the Af,-polarization measurement we 
already want to use just the semileptonic decays, so this alternative 6 tagging is actually 
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a high-efficiency option. We focus on 6 —^ rather than on 6 —?■ e decays because of the 
cleaner identification of muons in a typical hadron-collider detector. 

We estimate the performance of the “soft muon” b tagging using the CMS public 
note [65]. Ref. [65] gives the efficiency of selecting b jets versus the rejection rate for non- 
b jets, separately for selections based on and for the impact parameter. Requesting 
a large impact parameter is somewhat more effective than requesting large against 
light-flavored jets. However, this is largely due to the contribution from the 6 —>■ c — 
decay chain, which in our case is not part of the signal. As an example working point 
we therefore choose a p^^'-based selection that gives a survival probability of eudsg — 0.3% 
for jets initiated by u, d, s quarks or by gluons, Cc 2.5% for jets initiated by c quarks, 
and €b = 8% for true 6-flavored jets. The value for eb is inclusive, encompassing b quarks 
decaying directly into a muon, 6 —)> p, decaying indirectly, 6 —)• c —)• /x, or with no muon at 
all in the decay chain. For jets whose initial 6-quark decay chain does contain a prompt 
muon (6 —)• c/u/) the tagging efficiency is approximately 70%. 

4.2.2 Xc reconstruction 

The inclusive Xc state in Ab —?■ Xc^iv usually contains a Ac baryon, which often decays 
into a A. We examine three Ab selections in order of decreasing statistics, but increasing 
purity, 

• Inclusive selection: requiring only the presence of a soft muon inside a jet, 

• Semi-inclusive selection: requiring in addition the presence of a A —)• pvr” candidate, 

• Exclusive selection: completely reconstructing a Ac candidate in a few clean modes 
with only charged particles in the final state. 

The studies of A;, in QCD events [17-19, 28] use A;, —)■ J/'i/; (—)• p+/x“) A (—)• pvr“). Us¬ 
ing this decay would provide a cleaner sample than the three approaches described above, 
but it has a very small branching ratio of ~ 3.2 x 10“^. This requires large statistics, mak¬ 
ing it prohibitive to use in new-physics processes. Another clean decay used by LHCb [66], 
Ab —>■ n~)pK~, likely also suffers from a small branching ratio (not yet re¬ 

ported). However, these decays can become useful in the future for cross-checking and 
refining the information obtained from SM calibration samples like tt using the semilep- 
tonic decays on which we focus here. 

To measure the A;, polarization it is necessary to reconstruct the neutrino and A;, 
momenta. This is equivalent to knowing the XcP and neutrino momenta. In this subsection 
we explain how the XcP momentum is determined in each of the above selections. In the 
next subsection we use this information to also obtain the neutrino momentum. 

Inclusive selection In this approach we only require the soft muon, which for the signal 
events originates from A^ decays. (There is also a small contribution from polarized 'Eb 
baryons, discussed in appendix B.) An important background is semileptonic B decays. 
Even though these decays are isotropic, their presence in the sample dilutes the observables 
sensitive to the A^ polarization, cf. eq. (4.6). In the inclusive selection the purity of the 
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sample is small, 0(/baryon), as the branching ratios of semileptonic A;, and B decays are very 
similar. On the positive side, the data set is very large. As an estimate for i3(Afe —)• Xc^iv) 
in our numerical estimates we shall use B{Ai, —)• Adiv + anything) ~ 10% [51], neglecting 
the small contribution from decays in which contains a D meson rather than the Ac 
baryon (see discussion below). 

In this inclusive approach the A;, four-momentum can be determined only approxi¬ 
mately, and on a statistical basis. For b quarks produced at energies near the electroweak 
scale, the A;, carries on average only {z) ~ 70% of the 6-quark energy, with a broad distri¬ 
bution [67-70]. Approximating the z distribution with its average, we write 

~ {z)Eb . (4.8) 


To estimate the XcH energy, we first correct the measured jet energy, Fjetj by subtracting 
the energies of charged tracks originating from the primary vertex (assuming they are vr^) 
to obtain E'^^^. To get the Xc/r energy one would need to also subtract the energy of neutral 
particles from the primary vertex (mostly due to vr'^’s), Flneutrai) 


EXcbi- — '^jet A'neutral ■ 


(4.9) 


However, Fineutrai cannot be experimentally distinguished from neutral particles from the 
Afc decay. We thus make an approximation; the probability for a pion to be a is 1/3, 
so on average 

-Fneutral — - E}, . (4:-10) 

Using eq. (4.8) and E/^^ = Ex^fi + Ei, we express Ex^^ in terms of the corrected jet energy 
and the yet-unknown neutrino energy as 


Ex. 


cl^ 


^z)Ei^, - (1 - {z))E, 
2 {z) + 1 


3 (^) p/ 
2 {z) + l j"*' 


(4.11) 


In the last step we neglected the E^, term since it is typically an order of magnitude smaller 
than the first term. The same procedure works for the background decays B —)• Xc^iv- 
We also need to determine the momentum, Px^p-- While the muon is readily identifi¬ 
able and measurable, the momentum of Xc requires additional approximations. In cases 
where Xc contains just a (ground-state or excited) charmed hadron, the direction of Px^ 
can be taken as the direction of the track-based jet it produces and its magnitude can be 
determined from Ex^ assuming mx^ — It is not crucial to use the precise c-hadron 

mass since the parent 6-hadron mass is relatively large. If Xc contains additional charged 
hadrons, typically pions, observed as tracks originating from the 6-hadron decay vertex, 
their momenta can be included trivially, and their energies subtracted from Ex^ (assum¬ 
ing they are vr^) to obtain the charmed-hadron energy. More problematic are neutral 
hadrons, which contribute energy to the jet but do not leave tracks. One cannot deter¬ 
mine whether they come directly from the 6-hadron decay (in which case they need to be 
treated separately) or only from the subsequent c-hadron decay (in which case they are 
included in mx^ — Th® former case is problematic. It is not very common since 
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usually Xc is a single charmed hadron — B(Ai, —k Ac^v) ~ 0.7 B{Ah —)• Ac/w^ + anything), 
13{B^ ~ 0.8 ^(5° ^ + anything), 13{B^ ~ 0.8 B{B^ 

£)(*) 

fiv + anything) [51], and only a fraction of the remaining decays are expected to con¬ 
tain TT^’s. However, despite the small size of these contributions, misreconstruction of such 
events in the background can potentially contribute a large bias to the measured Afb, 
considering the low signal-to-background ratio of this inclusive selection. 

One possible handle for reducing the background from B ^ D decays is the shortness of 
the Ac lifetime relative to the H-meson lifetimes: ~ 5, Tj^± ~ 2.5, ~ 2. 

This is even more significant than the difference between D- and H-meson lifetimes, which 
is already being used by ATLAS as one of the handles for tagging c jets while rejecting h 
jets [71, 72]. For example, the loose operating point from ref. [71] provides 95% efficiency 
for c jets with a factor of 2.5 rejection of b jets. Perhaps an analogous technique could be 
used in our case for accepting Ac’s while rejecting a significant fraction of D mesons. While 
designing the relevant algorithms or estimating their expected performance is beyond the 
scope of this paper, we encourage further work along this direction and note that the 
shortness of the Ac lifetime has already been used for background reduction in a A;, study 
by DO [73]. 

It may also be possible to estimate the background contribution and subtract it. One 
could, for instance, use high-pT’ b jets from QCD events as a control sample. In this 
case the 6’s have no longitudinal polarization, so that the measured Afb will be entirely 
due to misreconstruction. One could further improve the accuracy of the background 
prediction using embedding: for the process of interest one would first select a sample of 
fully reconstructed 6-decay events, and then replace the 6 jet with a kinematically equivalent 
semileptonic 6 jet from the QCD sample (with the momentum determined from the rest 
of the QCD event). Here we do not pursue these ideas further but rather consider less 
inclusive selections that significantly suppress the background contributions, while keeping 
the overall statistical uncertainties comparable to that of the inclusive selection. 

Semi-inclusive selection The large background from semileptonic B decays can be 
reduced by requiring among the jet constituents both a soft muon and a A baryon. In the 
vast majority of A;, —?• decays we expect the Xc to contain a Ac.^ We can then use the 
decay chain Ac —^ A(—> P'^~) + X, with B{Ac —> A -|- X) ~ 0.35 and B{A —> pvr“) ~ 0.64. 
Requiring a reconstructed A —)• pvr” decay inside the b jet and originating from the vicinity 
of the displaced vertex will eliminate most of the H-meson background. Some H-meson 
contamination may still be present due to Kg —)• 7r'*'7r“ decays mimicking A —)• pvr”. This 
can be suppressed with a modest efficiency loss by requiring that the invariant mass of the 
two tracks, if assumed to be pions, is incompatible with the known Kg mass [28]. 

A reconstructed A in the jet can also be used for reducing the background from A;, —)• 
Ac + X with the Ac decaying semileptonically. While in principle the sign of the lepton 
eliminates this background, this requires knowing whether the jet originated from a 6 or a 6. 

^ Experimentally, for example, B{Kh —>■ D°p'K~) ~ 0.1 x B{Ab —>■ Aj"7r“) [51]. We expect A;, —>■ D^pl~v 
to be suppressed relative to At, —>■ A^l~v by a similar factor and thus A;, —>■ XcK to be dominated by 
Ab —>■ Ac^Py + X decays. 


- 19 



Decay mode 

Branching fraction 

A[!“ —7- pK~ 7r~^ 

6.7% 

A+ —)• Att^ —)• p7r'''7r“ 

0.9% 

A+ —)■ pKs —)• p7r'''7r“ 

1.1% 

A]!" —)• Att^tt+tt” —>■ p7r'*“7r'*“7r“7r“ 

2.2% 

A]!" —)• pKsTT~^TT~ —)• p7r'*“7r'*“7r“7r“ 

1.2% 


Table 2. Branching fractions of the main all-charged decays of Ac. For A+ —)■ pK~TT~^, we used 
the average from ref. [26] dominated by the recent Belle measurement [76] instead of the much less 
precise PDG value of (5.0 ± 1.3)% [51]. The ratio of the two values was used to rescale the other 
branching fractions from their PDG values, since they were measured relative to A+ —>■ pK~TT'^. 

Sometimes this information is available from the rest of the event, e.g., in a reconstructed 
tt sample from the sign of the lepton in a leptonically decaying top [74]. If not, one can 
use the relative signs of the lepton and the A decay products. 

In the numerical estimates we will assume an ca — 30% efficiency for A —)• p7r~ 
reconstruction. This is larger than the efficiency of 10-16% quoted by CMS in ref. [28] 
because we believe that quality cuts can be relaxed. The maximal achievable efficiency 
is limited by tracking efficiency, which is around 60%, considering the pair of tracks in 
A —)• pTT~ and integrating over the cr distribution of the A [75]. It should be noted though 
that the installation of new tracking detectors in ATLAS and CMS in the next years will 
likely significantly improve the reconstruction efficiency of long-lived resonances like the A 
baryon. 

In the semi-inclusive selection, reconstruction is approximate, performed using the 
same procedure as for the inclusive selection. 

Exclusive selection In this approach the hadronic system Xc is reconstructed very 
precisely by first reconstructing Ac from its decay products in one of the channels where all 
the products are charged, and then adding charged particles whose vertices are compatible 
with the reconstructed Ac origin. The strongest point of this approach is that one can 
obtain the Ac four-momentum without approximations. Xc is then known completely if 
the Ac is accompanied only by charged particles, and is known approximately if there are 
neutral particles like vr^. Moreover, purity against B mesons is expected to be high. All 
this comes at a moderate cost in statistics. Table 2 summarizes the most promising Ac 
decay modes. The dominant one has B{A'^ —?■ pK~7r~^) ~ 6.7% [26]. CDF have already 
used this channel for studying the A;, [27]. Second in size are the modes with an additional 
vertex from A or Kg decays, which have a total branching ratio of around 5.4% [26, 51]. DO 
have already used one of these channels (A)!' —)■ pKs —)• pvt+tt”) for studying the A;, [73]. 
It may be noted that restricting the analysis to just a fraction of the Ac decays does not 
invalidate the inclusiveness assumption in eq. (4.1) on which eq. (4.4) relies, as long as all 
Afc —)■ Ac -b Ai decays are included. 

The reconstruction efficiencies achievable for the decays in table 2, which involve be¬ 
tween three and five charged particles in the final state and in part of the cases an inter- 
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mediate resonance, should be estimated with a detailed detector simulation. We note that 
CMS has reported 33% efficiency for the three-prong decay —)■ J/V’(— >• ^~)K^ for 

> 30 GeV [77] and 20% efficiency for the four-prong decay —)■ )• 

K^K~) for 23 < < 50 GeV [78] in inclusive QGD production. For —)■ 

decays in VF -|- c production, GMS had about 11% efficiency for > 25 GeV [79] and 

ATLAS had 32% efficiency for p-^^ > 8 GeV [80]. In the following, we therefore assume 
that on average an efficiency of ~ 25% is achievable for the Ac reconstruction. 


4.2.3 “Soft neutrino” reconstruction 

Knowing the four-momentum together with the flight direction of A;, suffices to de¬ 
termine the soft-neutrino momentum up to a two-fold ambiguity [81] (see also ref. [82]). 
Experimentally, the A;, flight direction is the direction between the primary vertex and the 
secondary vertex associated with the soft muon. The neutrino momentum perpendicular 
(parallel) to the A;, flight direction, {Pu), is 

P,^ = -P±, P!=-a±Vb, (4.12) 


where 

— m? — 2P'j)P\\ ^ — rn? — 2Pj)^£'^ E‘^P‘j_ 

““ 2 ( 7^2 _ ^ 2 ) ’ 4 ( 7^2 _ ^ 2)2 7^2 _ ^2 ' 

Here, Pj_) -P||) E and m are the XcP system’s momenta perpendicular and parallel to the 
Afc flight direction, its energy, and its invariant mass, respectively. Eq. (4.12) gives two real 
solutions for pJJ if 6 > 0, and two complex solutions if 6 < 0. We propose to discard events 
with complex solutions since the backgrounds are more likely to have negative b values. 
The two real solutions can be treated on equal footing, as in refs. [81, 82], as both carry 
information on the neutrino momentum although with different resolution. However, we 
illustrate in section 4.3.2 how for the tt example the full-event information can be used to 
solve the ambiguity. 

The precision of the neutrino reconstruction is limited by the uncertainty on the di¬ 
rection between the primary and the secondary vertex. The angular uncertainty is 


5a ~ 


5x 

7A, CTA, ’ 


(4.14) 


where 5x is the uncertainty on the relative position of the two vertices and 7 Aj is the boost 
factor. It should be compared with the typical angle a of the neutrino momentum, which 
for 7A^ ^ 1 is, very roughly. 


This gives 


Pu ^ (-Pi^)rest _ ^ 

Pu Pu Px/^ 7Ai, 


5a 

a 


~ 0.2 


5x \ 
50 pm J 


ta„ 


1.45 X 10 


-12 


-1 


(4.15) 


(4.16) 
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independent of 7Aj,. Even though this uncertainty is non-negligible, it is not prohibitive. 

The method outlined here is applicable to any sample of b quarks. In cases where the 
rest of the event does not contain invisible particles, e.g., pp ^ Z ^ bb, pp —?• bb with the 
second b in the event identified as decaying hadronically, one can also use the measured 
^rp as input to reconstruction. 

4.3 Measurement in pp — )• ti events 

We now apply the general strategy for measuring the polarization to pp —?• tt events. We 
estimate its sensitivity, under several simplifying assumptions, for 100 fb“^ at 13 TeV LHC. 
Performing such an analysis in ATLAS or CMS would be very useful for calibrating the 
6-quark polarization measurement. Given the approximate universality of the polarization 
retention factor r^, see introduction and appendix C, the value extracted from the tt sample 
would be an important input when measuring the polarization of 6 quarks produced in new- 
physics processes. 

The analysis strategy that we propose consists of the following steps: selection of a 
ft-enriched sample by requiring an isolated lepton and at least four high-pT jets; recon¬ 
struction of a Afe candidate; global event interpretation in terms of jet-parton assignment 
and reconstruction of the neutrinos by the exploitation of kinematic constraints; measure¬ 
ment of the forward-backward asymmetry of the soft neutrino in an opportunely chosen 
rest frame. 

4.3.1 Event selection 

The best compromise between statistics and selection purity is achieved by targeting the 
final state with a single isolated electron or muon from IT-boson decay, for which the total 
branching ratio is approximately 30%. Final states with two isolated leptons give better 
selection purity but the branching ratio is six times smaller; an all-hadronic selection could 
achieve a reasonable selection purity only by imposing very tight kinematic thresholds. 
An additional benefit of the single-isolated-lepton sample is that one can veto the decay 
chain 6 —cX —)> iuX' using the relative sign of the isolated lepton from the W boson and 
the non-isolated lepton from the A;, (see section 4.2.1), in conjunction with global event 
interpretation (section 4.3.2). 

As an example, we take as baseline the same selection as in ref. [83], a tt analysis in 
the single-leptonic final state based on about 20 fb“^ of 8 TeV data, in which traditional 6 
tagging is not applied. This analysis requires exactly one isolated lepton with pT > 26 GeV 
and \7]\ < 2.1(2.4) in the muon (electron) channel, and at least four hadronic jets with 
Pt > 30 GeV and \r]\ < 2.4. In this way, 208(230) x 10^ events in the muon (electron) 
channel are selected, out of which 86(100) x 10^ are estimated from detailed simulation to 
be genuine tt events. Most of the background is composed of IV-l-jets events, with smaller 
contributions from multi-jet QCD production, Drell-Yan, and single-top processes. 

Going from 8 TeV to 13 TeV collisions, the tt cross section increases by a factor 3.3 [84]. 
If we assume similar selection efficiencies as for 8 TeV we expect 1.4 (1.65) million tt events 
in the muon (electron) channel for 100 fb“^ of integrated luminosity. The cross section for 
the main background, the inclusive W-boson production, increases by 1.9 (as calculated at 
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Selection 

Expected events 



Baseline 

3x10^ ti+0{ltf‘) bkg 



Soft-muon b tagging 

5x10^ ti+0{lt)^) bkg 



Signal events (t —>• 6 —)• Aft — )■ pi'Xc) 

Purity (example) 

XAfb/-A-fb 

Inclusive 

34400 

0(/baryon) (e.g., 7%) 

±7% 

Semi-inclusive 

2300 X (eA/30%) 

70% 

± 8 % 

Exclusive 

1040 X (eA,/25%) 

30% 

100% 

±19% 

±10% 


Table 3. Approximate number of expected ti events surviving different selections in the Af, po¬ 
larization analysis, for 100 fb“^ at 13 TeV. Baseline selection indicates the request of exactly one 
isolated lepton (electron or muon) and at least four jets, as in ref. [83]. ca is the efficiency of 
A —>■ p7r~ reconstruction, the efficiency of A^ reconstruction in the channels of table 2. Events 
originating from both b and b are included in all numbers. In the last column, the expected statisti¬ 
cal uncertainty on the soft-neutrino asymmetry for the different selections described in section 4.2.2 
is reported assuming the indicative purity in the third column and = 0 . 6 . 


NNLO with FEWZ 3.1 [85, 86 ]), but there are large uncertainties on the fraction of events 
with four associated jets above the px threshold. 

After the soft-muon selection of section 4.2.1 is applied to the events passing the base¬ 
line selection, we expect roughly 540 000 tt and 17000 single-top events (mostly tW) [87] 
to remain in the Run 2 dataset. Here, the yields for isolated-muon and isolated-electron 
channels have been summed. The rejection of non-top backgrounds depends on the poorly 
measured fraction of heavy-flavored jets associated with W, Z and 7 * production. Taking 
the associated jet multiplicity and heavy-flavor compositions of these samples predicted by 
MadGraph [ 88 ] with standard settings, and assuming that the multi-jet QCD background 
can be neglected, we expect less than 30 000 background events. 

The above estimates can be viewed as conservative. One can increase statistics by 
adding a “soft electron” b tagging. Moreover, the global event interpretation, outlined 
in the next subsection, can be used to further increase the signal-to-background ratio by 
selecting mass windows around the nominal masses of the reconstructed top-quark and 
IT-boson candidates. In the rest of the section we therefore simplify the discussion and 
ignore all non-top and single-top processes, focusing completely on the true tt events. The 
expected numbers of events are summarized in table 3, in which we also list the expected 
numbers of events after the three approaches to XcP reconstruction. 

4.3.2 Global event interpretation 

The Aft — 7 - Xcpv reconstruction procedure described in section 4.2.3 determines the soft- 
neutrino momentum, and correspondingly the Aft momentum, up to a two-fold ambiguity. 
This ambiguity can be resolved by checking which of the two hypotheses is more consistent 
with the kinematics of the full tt event, since the reconstructed 6 -quark momentum and 
the missing energy that would be attributed to the hard neutrino from t —)• 1 T 6 —)• ii^b, 
differ between the two solutions. 
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The global event interpretation is also useful for vetoing events in which the soft muon 
and the soft neutrino come from a 6 —)■ c —)• cascade. Such events can be rejected by 
requiring that this muon has the same (opposite) sign as the hard lepton coming from the 
opposite (same) reconstructed top quark in the event. This is important mostly in the 
inclusive approach to Xc reconstruction from section 4.2.2, where the charges of the Xc 
constituents are not measured. 

There exist various approaches to kinematic reconstruction of events with tops (e.g. 
refs. [89-95]). An important issue is that standard algorithms misreconstruct the ti event 
in a large fraction of the cases. For example, a radiation jet sometimes provides a better 
fit to one of the nominal ti products than the actual corresponding jet, especially when the 
latter is mismeasured or falls outside of acceptance. While extensive simulation would be 
necessary to determine which algorithms are best in our context and what their performance 
is, we would like to make several remarks. 

First, the reconstruction does not need to be fully correct for our purposes. In par¬ 
ticular, a correct reconstruction of just the top quark that produced the A;, suffices for 
resolving the soft-neutrino ambiguity and for vetoing events with wrong-sign leptons. It 
may even be beneficial in some cases to focus on reconstructing the relevant top rather than 
insist on reconstructing both. Also, even when the event is completely misreconstructed, 
the soft-neutrino solution will still be correct (accidentally) in roughly half of the cases. 

Second, we note that in the standard ti reconstruction approaches, the possibility that 
a significant fraction of the 5-quark momentum is carried by a neutrino is not taken into 
account. The prevalence of such events degrades the overall resolution of the reconstruction. 
Since we account for the soft neutrino explicitly, the reconstruction in our case will profit to 
some extent from this, usually ignored, additional information. The resulting impact on the 
performance of event interpretation depends on the applied algorithm and its estimation 
is beyond the scope of this paper. 

In the sensitivity estimates below, we will optimistically neglect the potential impacts 
of misreconstructed events. However, note that even if the ti reconstruction were com¬ 
pletely useless (which is an unreasonably pessimistic assumption), one could keep both 
soft-neutrino solutions and account for this ambiguity when interpreting the results, as 
was done in refs. [81, 82]. 

4.3.3 Expected sensitivity 

After resolving the ambiguity in the soft-neutrino momentum as outlined above, we apply 
the asymmetry analysis of section 4.2. In the last column of table 3 we collect estimates 
for the purely statistical component of AAfb/-^fb that follow from eqs. (4.6) and (4.7), 
assuming as an example = 0.6, c.f. eq. (3.25). These will also be the statistical uncer¬ 
tainties on the value of extracted from these measurements. 

We see that despite the different degrees of inclusiveness the three selections have 
comparable statistical uncertainties. Therefore, the fully inclusive selection is disfavored, 
considering the background reconstruction uncertainties discussed in section 4.2.2. The 
amount of background in the semi-inclusive approach is much more manageable, although 
the measurement would still be somewhat limited by the systematics related to the approx- 
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imations made in the reconstruction of the 4-momentum in the signal. The vertexing 
uncertainty described in eq. (4.16) is common to all the three approaches. Since many of 
the systematic uncertainties depend on experimental details that are difficult for us to sim¬ 
ulate using publicly available tools, and since the hrst measurement will likely be limited 
by statistics, the detailed study of systematic uncertainties is deemed outside the scope of 
this work. 

Overall, this looks like a promising measurement for Run 2 of the LHC. 

5 c-quark polarization measurement via A+ —)■ pK~7r~^ decays 

In principle, the semileptonic decays of Ac are similar to those of A;,. In this case it is the 
charged lepton rather than the neutrino that has approximately maximal spin-analyzing 
power. Unfortunately, the semileptonic channel seems impractical. First, its branching 
ratio is small, B{Ac —)• XpLv) ~ 3.1% — this estimate follows from rescaling B{D^ —)• 
Xfi^u) by the ratio of A^ and lifetimes. At the same time, semileptonic D decays, 
which constitute an intrinsic background, have much larger branching ratios, by factors 
of about 5 and 2 for and , respectively. This is different from the A;, case where 
the semileptonic branching ratios of B mesons and A^ are similar. Another difficulty is 
that, due to the relatively short lifetime, ~ 2.0 x 10“^^ s [51], there is a prohibitively 
large uncertainty on the Ac flight direction reconstructed as the direction between primary 
and secondary vertices, cf. eq. (4.16). Also the uncertainty due to additional neutral 
hadrons produced at the primary vertex is larger since they carry a larger fraction of the 
c-quark momentum than in the 6-quark case. Finally, backgrounds with prompt muons 
become more significant. The reason is the short Ac lifetime and the small rric', they make 
selection techniques that use impact parameter and relative muon pT, respectively, much 
less effective. 

A more promising decay mode is A+ —>■ pK~TT~^. Its branching ratio is relatively large, 
about 6.7% [26], while the H-meson background can be reduced significantly, without losing 
much signal, by restricting the invariant mass of the three candidate tracks to lie close to 
the Ac mass. For the angular distributions of each of the decay products we expect the 
same functional form as in eq. (4.2), but theoretical uncertainties on the hadronic matrix 
elements preclude precise predictions for the values of Oi for p, K~, and tt"*". These can 
be measured in the SM calibration sample. It may be noted that they can have different 
values for the different processes contributing to the pK~TT~^ hnal state, which include 
pK*(892)^, A(1232)++Ar“, A(I520)7r+, and non-resonant production [21]. Results from 
the NA32 experiment [20] indicate that aj^- is 0{1), as was conjectured in ref. [96], while 
Up and a^+ are small. 

5.1 Strategy for Ac-polarization measurement 

A way to tag a c jet for the purpose of the polarization measurement is to demand the pres¬ 
ence of a candidate A+ —)• pK~n~^ decay and consistency with a global event interpretation. 
As an example for the latter, we discuss in subsection 5.2 the polarization measurement of 
c quarks from W decays using a tt sample (see figure 4). 
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The identification of —)• pK~ 7r~^ decays in ATLAS and CMS is not trivial because 
the identities of charged hadrons are typically not determined by the detectors.^ We, 
therefore, propose the following strategy. Select three candidate tracks based on lifetime 
and vertexing criteria, i.e., requiring incompatibility with the hypothesis of tracks origi¬ 
nating from the primary vertex and compatibility with the hypothesis of coming from a 
common secondary vertex. The kaon candidate is the track whose charge is opposite to 
the other two. In some scenarios, the global event interpretation would tell us whether we 
expect a A^ or a A^, , and then a requirement on the absolute charges of the tracks can 
be added to reduce the background. Among the remaining two tracks, the one with the 
higher momentum (in the lab frame) is taken to be the proton candidate, and the other 
the pion candidate. This is almost always the correct choice for high-pT’ Ac’s because the 
proton is much heavier than the pion. In the small fraction of cases where this assignment 
is incorrect, the reconstructed Ac mass will typically fall outside the expected range, so the 
contamination will be minimal. After this identification procedure, the forward-backward 
asymmetry Afb of any of the three decay products {p,K,'k) in the Ac rest frame can be 
used for the polarization measurement. 

Since both the A^ —?■ pK~ Tr~^ branching fraction, table 2, and the Ac fragmentation 
fraction, eq. (2.3), are small, the intrinsic background from Ac decays to other final states 
(e.g., A^ pK~TT~^TT^, S'’'7r“7r'’', 7r^7r“7r'’“A) and D-meson decays (e.g., D~^ —?■ 

TT^ —)• TT~^K~7r~^7r~; Df —?■ K~^ K~'k^ , is a concern even after 

demanding the invariant mass of the p, K, vr candidates to be consistent with m\^. However, 
there are several effective handles for reducing many of these backgrounds: 

• In the signal decay, the kaon momentum in the lab frame is typically in-between the 
momenta of the pion and the proton, similarly to the discussion above. Demanding 
such a hierarchy reduces the background since in most of the background decays that 
contain three charged particles these particles are not p,K~,tt~^ so the negatively 
charged track does not necessarily tend to be intermediate in momentum. 

• Decays in which one of the three tracks, or an extra neutral particle, is a long-lived 
strange hadron, can be eliminated by vetoing on additional further-displaced vertices. 

• A veto on a fourth track consistent with the candidate Ac vertex can eliminate most 
of the backgrounds since the cannot decay to three charged particles. 

• The different lifetimes, r(A+, D+, D’*') ~ (2,4,5,10) x 10“^^ s, can be used for 

reducing all D-meson backgrounds. 

• Backgrounds from particular decays to three charged particles, such as —)• 7r'*'iL“7r+ 

and Df —)• K^K~Tr^, can be targeted by demanding that if the tracks are assigned 
the masses of these decay products, the resulting invariant mass should be inconsis¬ 
tent with that of the parent D meson. 

^Although particle-identification procedures based on specific energy loss or time of flight have been 
developed in both ATLAS and CMS [97-101], they show sufficient separation between protons and lighter 
hadrons only up to track momenta of D(GeV) at most. This is too small for the end-products of the decays 
of top quarks or new heavy resonances. 
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Figure 4. An example tt event that can be used for measuring the polarization of c quarks produced 
in W decays. 


• Mild cuts on the pt’s of the tracks and of the Ac candidate would be benehcial for 
reducing the background due to secondary hadrons produced in fragmentation (soft 
Ac’s and D mesons, as well as other sources of contaminating tracks). 

The combination of the above requirements will likely greatly suppress the backgrounds 
while reducing the signal by less than an order of magnitude. Yet, it is not clear whether the 
backgrounds will be negligible in the end. The precise amount of background is scenario- 
dependent. This is because the displaced-vertex properties, the ordering of the three tracks’ 
momenta in the lab frame for both the signal and the backgrounds, and the reconstructed- 
mass resolution, all depend on the charmed hadron momentum. The signal efficiency and 
purity will therefore depend on the kinematics of the process producing the c quarks. 
Estimating those for any particular process requires a detailed simulation and is beyond 
the scope of this work. In any case, since the Ac mass peak is narrow while the back¬ 
grounds are smooth, one can use a sideband for estimating and subtracting the bias that 
the backgrounds may be contributing to Afb- The background under the peak would still 
contribute statistical fluctuations. 

5.2 Measurement in pp —?■ tt events 

We now describe the Ac polarization measurement in pp —>• ti events, in which longitudinally- 
polarized charm quarks are produced via t —)■ W^h —)• csb, as illustrated in figure 4. We 
estimate the sensitivity for 100 fb“^ at 13 TeV under various simplifying assumptions. Per¬ 
forming such an analysis in ATLAS or CMS would be very useful for calibrating the c-quark 
polarization measurement. Such a calibration measurement is even more important than 
for Ab, because of possibly large Aqcd/?71c corrections to and r^, and the fact that the 
spin-analyzing powers of the A+ —)• pK~Tr^ decay are a priori unknown. 
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Selection 

Expected events 

Purity (example) 

AAfb/-^fb 

Baseline 

1.7x10® tt + C>(10®) bkg 



A+ —)■ pK~F~^ 

810 X (eAc/25%) 

20% 

100% 

26% 

11% 


Table 4. Approximate number of expected tt events surviving different selections in the Ac po¬ 
larization analysis, for 100 fb“^ at 13TeV. Baseline selection indicates the request of exactly one 
isolated lepton (electron or muon) and two jets passing standard 6-tagging selection out of at least 
four, as in ref. [102]. indicates the efficiency of Ac reconstruction in the A+ —>■ pK~TT~^ channel. 
Events originating from both c and c are included in all the numbers. The last column shows the 
expected statistical uncertainty on the forward-backward asymmetry of the Ac decay product with 
the highest spin-analyzing power assuming = 0.6, for two different assumptions regarding 
the achievable purity of the selection. 

The strategy that we propose here is similar to the A;, analysis from the previous 
section. It consists of selecting a tTenriched sample by requiring an isolated lepton and 
at least four high-p^ jets, reconstructing the event, and measuring the forward-backward 
asymmetry of the proton, kaon, or pion in the Ac rest frame. 

We start with a baseline selection of a single lepton and at least four jets similarly to 
section 4.3.1 and apply standard 6-tagging algorithms to remove most non-top background 
events. As an example we use the efficiencies from ref. [102], where the event selection 
contains a single isolated lepton (electron or muon) with px > 33 GeV and |r/| < 2.1, at 
least four jets with px > 30 GeV and \ri\ < 2.4, and exactly two of the four leading jets are 
required to pass a 6-tagging selection based on the combination of track-based lifetime and 
secondary vertices information. The 6-tagging working point corresponds to ej, = 70% [103]. 
With these selection criteria, 108 205 events survive in 20fb“^ at 8TeV with a composition 
of 94.3% tt, 3.4% single top (mostly tW), 1.9% lT-|-jets, and 0.4% Z-|-jets. We, therefore, 
expect roughly 1.7x10® tt events for an integrated luminosity of 100 fb~^ at 13TeV and 
we can neglect the non-top backgrounds. Event reconstruction can be performed similarly 
to section 4.3.2. Conventional 6-tagging algorithms can be used to assist the assignment of 
the jets. Ac candidates from the two jets interpreted as originating from 6 quarks should be 
vetoed. The expected number of signal events after reconstruction, using from eq. (2.3) 
and i3(A+ —^ pK~TT~^) ~ 6.7% [26], is shown in table 4. 

Let us estimate the expected sensitivity assuming that just one of the Ac decay products 
is being used in the polarization measurement, presumably the one with the largest spin¬ 
analyzing power cTj. Since it is likely that a* is close to 1 for the kaon [20, 96], and the 
possible values of rp are given by eq. (3.33), we will present estimates for atrx = 0.6. 
Considering the intrinsic backgrounds discussed in section 5.1, the signal efficiency eA^ and 
purity / cannot be determined without a detailed study. For the purpose of our estimates 
we assume = 25% as in section 4.2.2 and consider two possibilities for the purity /: 
100% and 20%. The resulting statistical uncertainty on the polarization measurement, 
^Afb/Afb, determined along the lines of eqs. (4.6) and (4.7), is shown in the last column 
of table 4. 

Overall, performing this measurement in Run 2 of the LHC seems feasible. 


- 28 - 








6 Isolating decays 


As discussed in detail in section 2, a large fraction of Af,’s are produced from the decays 




Afovr- 


±,o 


( 6 . 1 ) 


So far, we considered them part of the A;, sample. In principle, they can be distinguished 
from primary A^’s by observing a pion that together with the A^ reconstructs the ' mass. 
In practical implementations Q = m^Ahir) — m{Ah) — m-^ may be a better variable than 
the ' mass, because it reduces resolution effects from the A^ reconstruction. Vetoing 
the —?• Afovr^’*^ contributions would eliminate the leading depolarization effect, giving 

an even more direct correlation between A;, and 6-quark polarizations. In this section we 
discuss the prospects for identifying ' (and analogously Sc j decays at the LHC. 

An immediate difficulty is that the pion is very soft, m (*) — ttiaj ~ O.OdmAj,. In the 
semileptonic channels advocated in section 4 for the 6-polarization measurement, the A^ re- 
construction is not sufficiently precise for reconstructing S^ "^’s. This is due to the neutrino, 
whose reconstruction involves non-negligible uncertainties from the direction between the 
primary and secondary vertex, and due to the ambiguities surrounding neutral particles 
in the jet. Another difficulty is combinatorial background. The soft pion stems from the 
primary vertex, where additional pions and other hadrons are frequently produced as part 
of the jet in the 6-quark fragmentation process. In the case of a neutral pion, neutral 
hadrons produced in the A;, decay would contribute an additional ambiguity. It is thus 
likely that the optimal choice is to treat decayed "^’s as part of the A;, sample, as we 
have done throughout this paper. 

On the other hand, separate studies of ’ decays in the SM calibration samples, 
using well-reconstructed A^ decay channels where all the final-state particles are charged, 
could be very useful for better characterization of the polarization-loss mechanisms. The 
parameter A discussed in section 2 can be determined from the overall yield of these decays. 
wi can be determined either from the angular distribution of the pions (as discussed in 
ref. [11] and already attempted by DELPHI at LEP [40-42]) or from the Af, polarization. 
Overconstraining the system would even allow going beyond the dominant polarization-loss 
effects we consider in this paper. 

Using the formalism of section 3.2.2, we find that finite-width effects can be important 
in determining wi from the pion angular distributions. In the rest frame, relative to 
the direction of motion of the ^ in the lab, they are given by 


1 dT _ 1 9 

P dcosO 2 8 





( 6 . 2 ) 


In the narrow-width limit T (*) <C A, a = 0 for the T,^ (whose angular distribution is 
therefore insensitive to rci) and a = — 1 for the This case was emphasized in ref. [11] 
and assumed in the DELPHI measurement [40-42]. In the opposite limit, T (*) ^ A, 

a = —2. More generally, a depends on the reconstructed mass E of the SI ' as 


a{E) = —2 -|- 


8 (ms* -wrsj^ 


4[2{E - msJ2 + (E - ms*)^] + 3r2(E) ' 


(6.3) 
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Figure 5. Coefficient a from eq. (6.2) describing the angular dependence in Af,7r (thick 

blue curve) overlaid on top of the spectrum (black curve, arbitrary y scale). Dashed vertical 
fines show ±r (.) /2 ranges around the nominal masses, and horizontal fines indicate the values of 
a in the narrow-width limit. 

This is plotted in figure 5 for the masses and widths from table 1 and eq. (3.19). The value 
of a varies significantly within the and S), ^ peaks, and the average values within each 
peak may modestly deviate from the narrow-width-limit value, depending on the definition 
of the peak boundaries. More interestingly, the large widths provide an opportunity for a 
more precise measurement of wi. For example, one can focus the analysis on mass ranges 
with large |a| and/or avoid cancellations of sensitivity between mass ranges with positive 
and negative values of a as it happens in the case. 

Identifying S/ decays in ATLAS, CMS, and LHCb seems possible. One relevant A;, 
decay channel is 

Aft —)■ A/“7r“ , A/" —)• pA'^vr"'". (6.4) 

Using this channel, CDF has successfully studied s[*^’s at the Tevatron [32, 33], although 
not the quantities relevant in our context. LHCb has reconstructed Aft’s in this channel 
in ref. [30], although without reconstructing s[*^’s. In ref. [104] they studied —>■ Hftvr 

decays (which are analogous to —>■ Aftvr decays) using the channel —>■ H+vr”, H+ —>■ 

pK~ Tr~^. Another possibility is 

Ab^Jj-ipA, J/'i/—)• , A —)■ p7r“ . (6.5) 

This channel is already being used by ATLAS [18], CMS [19, 28] and LHCb [17] for Aft 
measurements. CMS has also studied —)• Hftvr decays using the similar channel H/" —)• 
J/'ijj —)• E~ —)• Avr”, A —)■ pTT~ [105]. Another possible channel, used by 
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LHCb in ref. [66], is 




( 6 . 6 ) 


For the decay chain in eq. (6.4), the spin-analyzing power is expected to be close to 
maximal [25, 106-111]. For the decay chain in eq. (6.5) there is disagreement between 
different theoretical approaches [25, 107, 110, 112-117], many predicting the analyzing 
power to be 0(0.1). The analyzing powers of the decay in eq. (6.6) are unknown. Not 
having a prediction for the spin-analyzing power is not a problem by itself since one can 
still extract wi from the polarization measurement by normalizing the result to a sample 
not enriched in '’s, or from the angular distribution of the pions from —)• A^tt as 

discussed above. 

For the c-quark polarization measurement, the idea of vetoing on Sc ' contributions is 
somewhat more promising. The pion is less soft, m (*) — ~ O.OOmA^, and the decay 

mode advocated in section 5, A+ —)• pK~7r~^, is fully reconstructible. Reconstruction of 
Sc*^’s in this channel has been performed by CDF in ref. [118]. 

The study of isolated Sc*^ samples is even more important than s[*^ since the informa¬ 
tion that can be obtained from inclusive Ac measurements is possibly limited. In particular, 
a direct measurement of may be problematic: the polarization in QCD events is sizeable 
only for momenta pc ~ ttIc, which is probably too close to Aqcd for factorization to be re¬ 
liable. Instead, one may prefer to use the theoretical prediction for (section 3.4), which 
relies on knowing A and wi. These two parameters can be obtained from measurements of 
the Sc yields and the angular distribution of the pion in S* —)• AcVr decays, respectively. 
The latter measurement has already been performed by CLEO [43], but it would be desir¬ 
able to verify its result, eq. (2.11), in view of the apparent discrepancies described at the 
end of section 2. Direct measurements of A and wi would also be useful for comparisons 
with theoretical models, considering that even in the longitudinal case the polarization 
measurement is only sensitive to the products aiV^ and the spin-analyzing powers ctj are 
unknown. It may be noted that since measurements of A and wi do not require a polarized 
sample, they can also be done in Belle, where high-precision studies of Sc ^’s have been 
reported recently [119], and in BaBar. 

7 Conclusions 

We pointed out that b and c-quark polarizations can be measured at the LHC, and de¬ 
signed general techniques that can be used for that purpose in ATLAS and CMS. The 
most interesting application would be characterization of new-physics processes producing 
such quarks. While new physics is yet to be discovered, we motivated a set of Standard 
Model analyses for ATLAS, CMS, LHCb, BaBar, and Belle that would help calibrate the 
polarization measurements. 

Our approach relies on the fact that A;, baryons partly preserve the initial 6-quark 
polarization. Since mb S> Aqcd; the processes that can change the polarization during 
hadronization are under good theoretical control. The dominant effect is due to ' de¬ 
caying to Afe and a pion [11]. While formally suppressed by Ijmb, the effect is numerically 
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0(1) for the values of and decay widths realized in nature. The depolarization 
effects can be parametrized by retention factors and ry for longitudinally and trans¬ 
versely polarized initial b quarks, respectively. Once ri and rx are measured in Standard 
Model calibration samples with known polarization, it will be possible to use them for 
studying the polarization of 6’s from possible new-physics processes. The same ideas apply 
to c quarks and the Ac baryons. 

Polarization measurements in Standard Model samples will also contribute to our 
understanding of QCD. As we discussed, there exist several different phenomenological 
approaches that give somewhat conflicting predictions for the non-perturbative QCD pa¬ 
rameters A and wi that determine ri and r^. Measurements of and vt in samples of 
quarks with a known initial polarization would thus be useful for assessing the ranges of 
validity of the various models. It would also be interesting to compare results obtained for 
bottom and charm quarks and examine to what extent the differences can be accounted 
for by higher-order effects in the HQET expansion. 

For a Ah polarization measurement, the semileptonic decay b —)• ciu seems particu¬ 
larly promising, with the neutrino being a perfect spin analyzer. For a Ac polarization 
measurement we suggest using A+ —)• pK~TT~^. 

We proposed to measure ri for b quarks using tt samples in ATLAS and CMS. After 
single-lepton tt baseline selection and identification of a potential A;, decay using soft- 
muon b tagging, the kinematics of the events is reconstructed. The 6-quark polarization is 
then probed by measuring the forward-backward asymmetry of the neutrino in the A;, rest 
frame. We examined several approaches, with varying degrees of purity, for dealing with 
the intrinsic background due to semileptonic B decays. In all of them, one can measure 
ri with about 10% precision using 100 fb“^ of data at the 13TeV LHC, considering only 
statistical uncertainties. While a full analysis of systematic uncertainties is beyond the 
scope of our work, we argued that at least in the high-purity approaches they are not 
prohibitively large. 

For measuring for c quarks, we again proposed to use single-lepton tt samples 
in ATLAS and CMS, where polarized c quarks are produced in W decays. Here, the 
calibration measurements will determine the products, of the c-quark polarization- 

retention factor and the spin-analyzing powers for each of the three decay products in 
A+ ^ pK- 'TT^. With 100 fb ^ of data, a precision of around 10%-30% is attainable. 

Finally, can be measured in the QCD production of 6 and c jets. As we discussed, 
and vt are different functions of several, currently unknown, QCD parameters. Therefore, 
measurements oi ri and are complementary. Reconstruction of A^ decays from which 
the polarization can be extracted, in inclusive QCD samples, was performed by LHCb [17], 
ATLAS [18] and CMS [19]. LHCb reconstructed also Ac decays [120]. We note that it will 
be useful for the polarization measurements in these samples to go beyond the constant- 
polarization ansatz assumed in [17-19] since the polarization is predicted to be a function 
of the parton-level kinematics of the event [16]. 

To reduce theoretical uncertainties it would be helpful to also have analyses that focus 
on Afe c’s produced from ^ decays. Besides studying the polarization of these samples, 
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(*) 

we argued that it would be useful to measure the yields (relative to the inclusive 

yields) and the angular distributions of the pion in ^ These analyses have to 

be done in fully reconstructible decay modes, where all the final-state particles are charged. 
Such studies can be performed by ATLAS, CMS, LHCb, and in the charm sector also by 
BaBar and Belle. 

Even though ri and are mostly universal, i.e., independent of the production mech¬ 
anism, they do have a weak dependence on the energy scale of the process. Their scale 
dependence is calculable by relating them to fragmentation functions. The required inputs 
can be acquired by measuring and at a fixed reconstructed 6-quark momentum but 
binned in the A;, momentum (and similarly in the c-quark case) once sufficient data are 
available. 

In short, the initial polarizations of 6 and c quarks are encoded in the polarizations of 
Aft and Ac baryons, respectively. The upcoming Run 2 of the LHC will allow measuring 
the universal retention factors with tt samples. 
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A More on Aft polarization for finite widths 

The results of section 3.2.2 were obtained by evaluating the integral in eq. (3.23) numer¬ 
ically. Here, we derive approximate analytic expressions by taking the energy-dependent 
factors PniE) and in eqs. (3.18) and (3.23) to be constant factors 

Vf and e-™/ 2 T 

in eq. (3.18). The energy-dependent widths r(E) in the propagators are replaced with a 
constant T, and the lower limit of integration in eq. (3.23) is set to —oo. 

We still need to specify which constant values to use for the m and T. For terms in 
eq. (3.23) that involve just the Eft or just the i.e., the non-interfering terms, it makes 
sense to take and ms*, respectively, for m and similarly Ps^ and Ps* for P (although 
the dependence on P does drop out after the integration in eq. (3.23)). For the interfering 
terms, on the other hand, it makes more sense to use some effective value rries between 
msj, and ms* and the corresponding Peff = r(mefr). 
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We then obtain 


V, = 


2R — I + 2 {1 + R) wi + 4i?efr (2 — wi) / + l) 


+ 


3(1 + 2ii) 

1 + i? — 2i?eff/ + l) 


rx = 


3(1 + 2R) 

1 + i? — 2i?eff/ (x^ + l) 
1 + 2R 


(2 — 3t(;i) sin^ 9p , 


wi — - j sin 6p cos 9p , 


(A.2) 

(A.3) 


where R has been defined in eq. (3.9) and similarly i?efr = e ^ 

A 


X = 


Teff 


(A.4) 


In the mi) —>■ oo limit, Sf, and have equal masses and widths, leading to R = Res = 1, 
X = 0, and thus Vz = ^, 'Px = 0, as expected. Since in reality x is 0(1), the deviation 
from the formal mi, —)• oo limit can be large. Even for R = i?eff = 1 the depolarization can 
still be 0(1). For instance, taking wi = 2/3, one has = (1 + llx^/27)/(l + x^). In the 
narrow-width limit, x —)> oo, eqs. (A.2)-(A.3) reduce to eqs. (3.13)-(3.14). 

Let us now substitute numerical values for Leff and i?efr in eqs. (A.2)-(A.3). For 
example, values corresponding to mes = give 

Pf'~ 0.17 + 0.41t(;i, Pj ~ 0.59 - 0.217X1, (A.5) 


while for mes = ms* 

P^^ ~ 0.28 + 0.367X1, Pj~ 0.64 - 0.187X1. (A.6) 


These numbers are close to the exact results in eq. (3.24), which lie between the two cases 
for meff. 

Another point we would like to make is that the physics of spin rotation, which we 
have been describing in momentum space, can also be described as oscillations in time 
between spin eigenstates, similar to K^~K oscillations, for example. The b spin in our 
case is the analog of strangeness, while the 'Eb and are the analogs of Kl and Ks- With 
the approximations made in this appendix, it is possible to interpret the physics in this 
way if we assume a common width T for the Eb and E^ and ignore the small effect of the 
thermal factor. With the Fourier transform 


1 

E — m + iT/2 


oc 


dt e 


iEt —imt—Vtl^ 


(A.7) 


for each of the propagators, we obtain 



dE\E) (P| oc 





(A.8) 


where 


|T™,(i)>oc ^P^.^(0p)^(i,M|4,+4; l,m)\Eb{M)) 


M 


+ e 


— 2777 .^* t, 


Y,Rm'm{Op)Y,{lM\l,+l- l,m)\El{M)) (A.9) 
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describes the time evolution (oscillations) of the state initially given by eq. (3.4). The 
time-dependent prefactor in eq. (A.8) describes the fraction of particles that decay at 
time t. 


B Sfc polarization 

In the main text we consider the most common hadronization of b into baryons, which is 
that b hadronizes with u and/or d quarks. However, in roughly 15% of the cases, one of 
the light quarks is s producing baryons. These are isospin doublets with spin 

configurations equal to the ones of A;,, respectively. The polarization formalism of 

section 3 thus applies also in this case. Polarized H^’s, produced directly, as well as from 
and decays, can be used to improve the statistics of the b polarization measurement. 

The mass splitting between and ~ 20 MeV, is much larger than their decay 

widths, Th* — 1.6 MeV and r=/ < 0.08 MeV [104, 105, 121]. The Ej, depolarization due to 
and decays can therefore be described in the narrow-width limit, eqs. (3.15)-(3.16). 
The statistical hadronization model gives in this case ~ 1.2 and R^i, — 0.91, consistent 
with partial information on the relative production rates [104]. The polarizations in the 
longitudinal and transverse cases are 

(pf)^^ ~0.10-F0.45t(;i, ~ 0.55 -0.23^1, (B.l) 

giving the total polarization retention fractions (after including direct Ef, production) 

~ 0.51, 0.67, 0.75, (rr)=^ ^ 0.75 , 0.67, 0.63 , (B.2) 

for wi = 0, 2/3, 1, respectively. These values of and vt are similar to those character¬ 
izing Ab, eq. (3.25). 

If the semileptonic decays of the b are used for the polarization measurement, the 
possibility discussed in section 4, one might consider performing an inclusive measurement 
including both A^ and Ef, contributions. The Ef, semileptonic branching ratios are expected 
to be similar to those of the Af, (for one of the dominant decays, Ef, —)• Ecii^ and Af, —>■ Acl'z/, 
see ref. [122] and references therein). Assuming equal branching ratios, and also that the 
system has the same value of wi as the system, the weighted averages that 

the inclusive measurement would be sensitive to are 

~ 0.45 , 0.64 , 0.73 , 0.73 , 0.64 , 0.59 , (B.3) 

for rci = 0, 2/3, 1, respectively. 


C Fragmentation functions for Ab and Ac 


In this appendix we express the polarization retention factors and vt defined in eq. (1.1) 
in terms of the 6 —)• Ab fragmentation functions. In particular, we want to show that 


riiz) 


Gijz) 

Di{z)' 


rxiz) 


Hi{z) 
Di{z) ’ 


(C.l) 
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where z is the fraction of the initial 6-quark momentum carried by the A;,, and the frag¬ 
mentation functions Gi{z),Hi{z)^Di{z) are defined below (see also ref. [123]). These can 
then be used to compute how and tt vary with the scale of the hard process. The same 
formalism applies to c —?■ Ac fragmentation functions. 

The cross section for a hadron h with transverse momentum pT and spin state Sh is 
given by, see e.g., ref. [31], 


da {pp —> h{Sh) -I- X -I- ...) 

dpT 


dpTdz Tr 


da {pp ^ q + ...) 

dpT 




iSh,z) 


5{pt - zpt), 


(C.2) 

where dd/dpT is the differential cross section for production of the hard parton q, without 
fragmentation, in the process of interest, X denotes the additional particles produced in 
the fragmentation of q, and ellipses denote all the other final-state particles. We have 
suppressed the dependence on the factorization scale p of both dd/dpT and A^, the latter 
containing the fragmentation functions. The fragmentation functions are universal, inde¬ 
pendent of the hard process. The trace in eq. (C.2) contracts the Dirac indices of the 
outgoing q in dd/dpx with those of A^. We are interested in the case q = b, h = Af, and, 
separately, also in the case q = c, h = Ac. 

In the ultra-relativistic limit, for a quark q with momentum k hadronizing to a spin-1/2 
hadron h with mass M/j, momentum and spin Sh, the relevant fragmentation functions 
Di{z), Gi{z), and Hi{z) are given by (see, e.g., refs. [38, 123-126]) 




^- ^^L^P^<{Q\q{i)\X-Ph,Sh){X-Ph,Shmm . (C.3) 


2 ( 27 r) 


«-=0 


where 


Xq{Sh, = \ {Di{z) i- + Gi{z) Ah 75 yl- + Hi{z) ia^^'y^n’tShj./Mh^ , (C.4) 


where we use light-cone coordinates with = (1, 0, 0, ±1) and take re_ to be aligned with 
k. The light-cone components of a four-vector are = a - nip. We also use z = P^/k~ as 
the light-cone fraction of the quark momentum carried by the hadron h and the sum is over 
all the hadronic states X that accompany h in the jet. is the spin vector describing a 
pure spin-1/2 state, which in the rest frame of the hadron is just Sj{ = Mh (0, Sh) (see, e.g., 
ref. [127]). It has been expressed above in terms of the light-cone helicity Xh = /P^ 

and the transverse components = 5^ — Sj^'nAL/2 + {Sf^ M‘^ / P^"^) n^l2. 

The fragmentation function Di{z) describes the probability for a certain hadron to be 
produced from a given quark. The total fragmentation fraction, like the ones quoted in 
eqs. (2.2) and (2.3), is then given by 


f,^h= I'dzDl^iz). (C.5) 

Jo 

The fragmentation functions Gi{z) and Hi{z) encode, in addition, the polarization of the 
hadron when produced from a quark with spin pointing in longitudinal and transverse 
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direction, respectively [123, 128]. We have suppressed the h, q indices on the fragmentation 
functions in eq. (C.4) and the fact that they depend on the factorization scale /i. 

For heavy quarks some control on the fragmentation functions can be achieved using 
HQET, see e.g. refs. [129-132]. This is relevant for the polarization of as the main 
depolarization effect indeed originates from the finite quark mass and can thus be described 
in HQET. In the exact —)• oo limit, the A;, spin is completely aligned with the spin of 

the b quark. Therefore, the product of matrix elements in eq. (C.3) has the same Lorentz 
structure as the outer product of the two Aquark Dirac spinors 

\ 1+i4/t 5'f, • eg ^ ^ bT \ 

urn = (1 + j = (1 - ^ ® ^^1 ■ ( C -®* 


Here, and Cg^ coincide with the hadron velocity four-vector, /Mh = {Eh, 0, 0,ph)/Mh, 

(where in our conventions ph < 0 ) and its longitudinal polarization vector, Cg = {ph, 0 , 0 , Eh)/Mh, 
respectively. They satisfy = 1, e| = —1, u • €3 = 0, and v ■ Sh = 0- In this formal limit 
the fragmentation function of a heavy b quark to h = Ab reads 


^iiSh,z) 


with 


2Mh 1+]J 
Eh-Ph 2 




Di{z) = Gi{z) = Hi{z) (X 5(1 — z). 


(C.7) 

(C. 8 ) 


We see that in the heavy-quark limit, the A;, fragmentation functions at /i < nib ai'e given 
by a single function. Eqs. (C.4) and (C.7) coincide in the ultra-relativistic limit in which 
yh = n^Eh/Mh -I- • • •, Cg = — n^Mh/2Eh -!-•••. Apart from RG effects discussed be¬ 

low, measurements with highly energetic jets thus probe deviations from eq. (C. 8 ), which 
are precisely the finite-mfe effects calculable in HQET. A perturbative treatment of heavy- 
quark fragmentation is possible, if the fragmentation function is summed over all possible 
final states [129]; for the nonperturbative endpoint region z ~ 1 see ref. [131]. We restrict 
ourselves to the exclusive case of fragmenting to one heavy hadron and make no assump¬ 
tions about the form of the fragmentation functions. Unpolarized fragmentation functions, 
Di{z), have been measured for inclusive samples of b hadrons at LEP [67-69] and SLD [70], 
and for the Ac by CLEO [133], Belle [134] and BaBar [135]. See refs. [131, 136] for theo¬ 
retical interpretations of such measurements. No measurements of polarized fragmentation 
functions are available yet. 

As argued in the main part of the paper, when departing from the heavy-quark limit 
the dominant effect of depolarization is due to the hadronization of the b quark into not only 
Ab but also into T,b and baryons. The T,/ "^’s decay to A^ via strong interactions, albeit 
with a phase-space suppressed decay width that parametrically enhances the depolarization 
effect. We have parameterized the relative production probabilities of the A;, and S), ' states 
using the nonperturbative parameters wi and A from eq. (2.7) and in the narrow-width 
limit R from eq. (3.9), to compute the polarization retention factors and tt- They are 
directly related to the fragmentation functions 


riiz) 


Gijz) 

Di{z)' 


rxiz) 


Hi{z) 
Di{z) ’ 


(C.9) 
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which are the two relations already advertised in eq. (C.l). The two relations can be easily 
understood from eq. (C.4) or (C.7). For instance, the longitudinal spin projector will select 
the Di ± Gi combination for positively (negatively) longitudinally polarized baryon, while 
Hi is similarly related to transverse baryon spin. 

For example, suppose that a Pl (Pr) projector in the hard kernel acts on the outgoing 
b quark, so that a left-handed (right-handed) b quark is produced. If we measure the spin of 
the Af, along its direction of motion, this projector gets multiplied by the linear combination 
Di =F 75 G 1 of the fragmentation functions in eq. (C.4). In the case that Di = Gi — 
like in the heavy-quark limit — the fragmentation function itself is proportional to the 
same projector; thus a fully longitudinally polarized A;, with negative (positive) helicity is 
produced after the fragmentation, i.e. vrIz) = 1, compatible with eq. (C.9). Oppositely, if 
Gi{z) = 0 and Hi{z) = 0, the probability of producing a Ab of specific spin would be 
independent of the underlying spin of the b quark; this is possible only if the polarization 
is lost completely, i.e. ri^Tiz) = 0. 

Now we would like to explain several points using the example of A^ production in 
e'’'e“ collisions at a specific center-of-mass energy, Scm- The cross section for producing 
a Aft with spin Sr is given by the usual convolution of hard kernels with fragmentation 
functions [131] 

(e^e~ Ab -F X) = f ^TT[Hi{Ecm, x, {Sh, z/x, ^)] . (C.IO) 

-i' ^ 


Here, Hi are the perturbatively calculable hard kernels, and x the fractions of the total 
available energy carried by the Ab and the initial parton i, respectively, and the sum is over 
all final-state partons. In principle one may worry about subleading corrections due to Ab 
fragmenting from an initial gluon or light quark. These corrections are process dependent 
and are practically negligible for the processes we are interested in. As an example consider 
Ab production on the Z pole, which is dominated by the e'’'e“ Z ^ bb partonic process. 
The longitudinal polarization retention fraction is given by 


Gijyjz) + Gf^i{z){al - (Jr)/{( rj - ^ r ) 

+ Ei i^) (4 + 4) / (4 + 4) ’ 


(C.ll) 


where the 4 h partonic cross sections for the left-(right-)handed b quark, and 

similarly for the other partons i = u,u, d, d,.... Here we see a small violation of universality 
not due to the scale dependence but due to the sum over light quarks and anti-quarks in 
eq. (C.ll). These non-universal contributions are suppressed by Q;^(mb) as they require 
fragmentation of a light quark (or antiquark) to a heavy-quark baryon, Dy\{z), Gf^iz), 
and are thus small (such perturbative fragmentation was calculated in ref. [137]). For Z 
decays, the fraction of events containing g ^ bb (or g ^ cc) is only about 0.3% (3%), 
as determined both theoretically and experimentally (see table 17.2 of ref. [51]). These 
numbers still need to be multiplied by the relative branching fraction of the total qq vs. bb. 

The next point we would like to make is that in the main text, vr and vt describe the 
average properties of the full sample, which includes baryons with all the possible values 
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of z. Thus, the retention factors are, for fixed center of mass as in the above example of Z 
pole or for t decays, given by 


f}dzGi(z) [} dz Hi(z) 

rL = ^ - rT = ^, -(C.12) 

fgdzDi(z) fgdzDi(z) 

Just like the fragmentation functions, they are independent of the production process, 
except for a logarithmic dependence on the hard scale. In the example of eq. (C.IO) 
integration over all possible z for i = b leads to 


= Tr / dxnb{Ecra,x,fi) 




z,d) 


(C.13) 


This demonstrates explicitly that as long as we are only interested in the polarization from 
Aft’s of all z’s only the inclusive retention factors are needed. 

The polarization retention factors and rx are universal, up to the logarithmic run¬ 
ning of the fragmentation functions with the characteristic energy scale of the process. 
Therefore the universality violations will be small if they are used for new physics mea¬ 
surements at scales not too different from the scale at which or equivalently Di, 

Gi, and Hi, are first extracted (e.g., in top decays as we propose in this paper). The 
fragmentation functions evolution is governed by perturbative splitting functions, similarly 
to the evolution of parton distribution functions (see, e.g., refs. [138, 139]). The resulting 
universality violations can be estimated for instance using the model for the fragmen¬ 
tation functions in ref. [38], with the LO RG running calculated in ref. [140] (see also 
refs. [123, 139]), in which the violation in rL,T universality is seen to be relatively mild. 
Taking the results of ref. [38] at face value the rx is found to change by 0(15%) due to the 
RG running between 5 GeV and 45 GeV for Aft (and by 0(10%) for Ac due to running from 
2.2 GeV to 45 GeV), while the change in rx is 0(5%) (0(15%) for Ac). We stress that 
these estimates apply only to the model of fragmentation functions as obtained in ref. [38], 
and could differ for the measured (in the future) shapes of fragmentation functions. 

Once sufficiently precise measurements of rx,xiz) are available, it will be possible to 
extract the fragmentation functions from them (using also information on unpolarized b- 
hadron production). The Aft production cross section and polarization retention in new 
physics models can then be calculated using factorization expressions as in eqs. (G.2) 
and (C.IO) after evolving the fragmentation functions to the relevant scale. 

Finally, we would like to comment on the experimental analyses [23, 27-30, 141], 
mentioned in section 2, which measured the px dependence of the ratios of fragmentation 
fractions for different b hadrons. Using eq. (C.2), restricting to q = b, the differential cross 
section for the production of an unpolarized b hadron is given by 


da^ 

dp!j^ 



dz dd^ 

Z dp^rp 



(C.14) 


The experiments report px dependences of the fragmentation fraction ratios, which are 
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thus given by 


dz da^ 


fb^hi (p: 


2 ; dpn 


Pt =Pt/^ 




fb^h2{PT) 


dz da^ 


dpi 


(C.15) 


Ou(^) 


Pt — Pt! ^ 


1,6 


We note that the dependence on the details of the hard process does not cancel out, as long 
as is not proportional therefore these ratios are not universal quantities. 

To extract the fragmentation functions, D^^z), and the fragmentation fractions, eq. (C.5), 
it would be useful to measure cross sections differential in two variables, in bins of both the 
reconstructed 6-quark pT and the reconstructed 6-quark momentum fraction carried by the 
Afc. The reconstructed 6-quark momentum is obtained by adding to the 6-jet momentum 
the momenta of neutrinos. 
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